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Interacting Self Avoiding Walk (ISAW) on a lattice is a simple model to study the 
Coil to Globule transition of linear homopolymers. The temperature at which the 
transition takes place is called the theta temperature. The value of theta temperature 
depends on the chosen lattices. The value of theta temperature for ISAW on FCC 
lattice is reported as 6.46 in a simulation study [J. Chem. Phys. 135 , 204903 (2011)]. 
Latter it has been reported as 7.614 [J. Chem. Phys. 138 , 024902 (2013)]. Simu¬ 
lations in these two studies involved long chains with over 2000 monomers. In this 
paper, we present a method that gives a reliable estimate of the theta temperature of 
ISAW on FCC lattice using chains of shorter walk lengths (less than 100 monomers). 
We compute density states of ISAW on FCC lattice by employing growth walk al¬ 
gorithm and then we use the recently introduced pseudo order parameter method 
to estimate the theta temperature. The value obtained from our method agrees 
very well with latter reported result of 7.614. In order to corroborate our result we 
also estimate the theta temperature using the celebrated partition function zeroes 
method. 
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I. INTRODUCTION 


Conformational transition of a linear homo-polymer chain from its coil state to globule 
state by a change in physical (temperature) or chemical (solvent concentration) conditions in 
dilute solution has been a widely studied phenomenal^. This transition is due to competition 
between excluded volume effect and the monomer - monomer interactions. These interac¬ 
tions get balanced at a particular temperature called theta temperature or at a particular 
solvent called theta solvent for which the linear polymer behaves as an ideal chairt^. The 
ideal chain behaviour has been studied by a random walk model of the polymer whose en¬ 
semble average of radius of gyration Rg scales with length of the chain N as < >~ 

where u is known as the Flory’s exponent which is equal to 1/2 for the ideal chain^^. 

The excluded volume interaction is included in random walk model through a constraint 
that segments of walk do not intersect with each other. Such a model is called as Self 
Avoiding Walk which is useful to study polymer in a good solvent condition or at high 
temperature^. Interaction between the monomers can be an attractive interaction among 
non bonded nearest neighbours. This is called as Interacting Self Avoiding Walk. The 
continuous^iSi^ lattice^‘i^‘^‘22i‘^ and off latticeii*^^ versions of this model have been studied 
extensively. 

Interacting Self Avoiding Walk (ISAW) on a lattice is a simple model to study the coil to 
globule transition of linear polymers. The transition temperature, known as theta temper¬ 
ature depends on the chosen lattice^. Since the system behaves as an ideal chain at theta 
temperature, it can be estimated based on the analysis of the value of Flory exponent u eval¬ 
uated from the scaling of ensemble average of the radius of gyration or end to end distance 
with walk length N. That is the temperature at which u = 1/2 has been identihed as the 
theta temperature!^. Several other methods can also be used to obtain theta temperature 
notably analysis of partition function zeroes!^ and inflection points^*^ of microcanonical 
entropy. In a recent study we have introduced a pseudo order parameter method to identify 
the theta temperature^^. 

Theta temperature of ISAW has been obtained for different types two and three dimen¬ 
sional lattices^'^*^'!^'!^'!^*^*^*^ . In particular, recent Monte Carlo simulation of ISAW on 
Face Centered Cubic (FCC)lattice reports an estimate of theta temperature as 6.46!^. The 
theta temperature reported in another work from same group is 7.614!^. These studies 
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employed Metropolis Monte Carlo method to compute Flory exponent as a function of tem¬ 
perature. The theta temperature is estimated as the value for which the Flory exponent 
is 1/2. Two different values of theta temperature reported on the same lattice initiated 
us to study ISAW on FCC lattice systematically. We compute density of states of ISAW 
on FCC lattice for the hrst time. We estimate theta temperature of ISAW on FCC lattice 
using recently proposed pseudo order parameter method and well known analysis of parti¬ 
tion function zeros on complex temperature plane. We hnd that a reliable value for theta 
temperature can be estimated from pseudo order parameter method using even ISAW of 
short walk lengths. 

The paper is organized as follows. Computation of density of states and microcanonical 
entropy of ISAW on FCC lattice for different walk lengths using the concept of atmosphere 
is presented in section II. Estimation of theta temperature using pseudo order parameter 
method^S as well as analysis of partition function zerosI^^^'I^i^Si*^ is discussed in section III. 
Conclusions are presented in section IV. 


II. DENSITY OF STATES ESTIMATES OF ISAWS ON FCC LATTICE 

In our simulation we employed growth walk algorithm^'^'^ to generate the self avoiding 
walk on a FCC lattice. There are many growth walk algorithms that can be used to compute 
density of states of ISAWs^^*^. In particular, the algorithm we have employed to generate 
SAWs are called Rosenbluth-Rosenbluth (RR) walk or Kinetic growth walk—. In RR algo¬ 
rithm, the walk start at arbitrary chosen origin on a lattice and search for the unoccupied 
nearest neighbors sites. The walk will randomly move to any one of the available unoccupied 
nearest neighbour site with equal probability. This process continues until the walk reaches 
a given walk length N. In walks, if there is no unoccupied nearest neighbour sites available 
for a walk to proceed further to reach given N, such walks are called trapped walks. These 
walks can be discarded and a fresh walk can be started from the origin. All the SAWs 
generated in this algorithm are rooted at the origin of the lattice. By assiging an energy e 
to all non bonded nearest neighbors (nbNN) of succesfully generated self avoiding walk of 
length N one can obtain ISAWs. The energy of a conhguration with ‘m' nbNN contacts 
are E = me. Without loss of generality one can take e = — 1 for attractive interaction. An 
ensemble of such walks are generated. 
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Statistical mechanics of a given system can described, once we know the Density of states 
(DoS). Counting Self Avoiding Walk conformations corresponding to a particular energy 
will in turn give the DoS for that particular energy. We used the concept of atmosphere in 
our growth algorithm and counted the DoS of ISAWs^^— . For a particular conformation C, 
let there be a, number of unoccupied nearest neighbors sited available at step, then a* 
is the local atmosphere for that step. Global atmosphere or simply the ’atmosphere’ A for 
the successfully generated N step walk is dehned as^*^, 

N 

A^(c) = n«*(c). (1) 

i=l 

The atmosphere for trapped walk is always zero^. Since FCC lattice has 12 nearest neigh¬ 
bors, highest number of nearest neighbors z in any other lattice in three dimension, the 
trapped conformations are considerably less in the FCC lattice. Let there are m{C) number 
of nbNN contacts in a particular conformation C, energy of that particular conformation 
E{C) is me or —m. Let denotes the Density of states (DoS) ISAWs of length N having 
energy -m, the Monte Carlo estimate of DoS can be obtained from the atmosphere as^*^ 

1 

^ E A^{C,) (2) 

1=1 

Here M is the total number of walks including trapped conhgurations. Mm is the total 
number of walks with energy —m. The monte carlo estimate of microcanonical entropy can 
be found from the estimated density of states as, 

S^{E = -m) = log{CfJ (3) 

Figure [1] presents the computed microcanonical entropy for ISAWs of different walk lengths. 
We have computed DoS up to walk length 300. In a small cluster, for 300 step walk, the 
computation took more than a month to complete generation of 10® walks. Hence we 
limited our further analysis for ISAW of walk lengths less than 100. We show that the 
theta temperature can be reliably estimated using recently proposed pseudo order parameter 
method using ISAW of short walk lengths. We observe however higher walk lengths up to 
300 are required for a reliable estimate of theta temperature using the analysis of partition 
function zeros. 
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FIG. 1. Microcanonical entropy of ISAWs on FCC lattice for different walk lengths 

III. ESTIMATION OF THETA TEMPERATURE OF ISAWS ON FCC 
LATTICE 


At a given inverse temperature {3 = where is the Boltzmann constant (which is 
taken as unity), the Canonical Partition function can be calculated from the DoS estimates 
as 

™maa; 

Cn{I3) = 51 C'Ar,mexp[-me/?], (4) 

m=0 

where is the maximum number of contacts that an N step walk can have. The canonical 
probability distribution is defined as^, 


PN{/3,m) 


Cat,^ exp [-me/3] 
CN{f3) 


(5) 


It may be noted that in the present work, e = — 1. It is well established that the 
coil-globule transition in polymers are second order transitions^^I^iS^. Canonical probability 
distribution as a function of number of contacts is shown in fig. [2] for two values of inverse 
temperatures /3 on either side of inverse theta {/3e) temperature. For second order transitions 
canonical probability distribution exhibits a single peak with the location of maximum of 
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FIG. 2. Canonical probability distribution at two different /I below and above /?© for ISAW of 
walk length 40. 

the peak varying continuously with (3 across /9©. In contrast, the canonical probability 
distribution shows a double hump for systems exhibiting hrst order phase transition and 
location of maximum of the distribution discontinuously jump from one peak to the other 
at transition temperature. This jump can be used to estimate the transition temperature in 
hrst order transitions. 


A. Pseudo order parameter method 

It has been observed recently that the canonical probability distribution for systems 
exhibiting second order transition shows a single peak with an asymetric long tail^S. The 
long tail jumps from one side to other side across the transition. Based on this observation 
a pseudo order parameter is dehned to obtain a binary signal of theta point transition on 
two dimensional square lattice^^. For ISAWs on FCC lattice, it can be seen from £g. [2] that 
long tail is on higher contact side for (3 = 0.108 < /S© and it is on the low contact side 
13 = 0.209 > (3e. 
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FIG. 3. Pseudo order parameter versus inverse temperature for different walk lengths 


The pseudo order parameter defined based on the shifting of the tail of canonical proba¬ 
bility distribution function across the theta temperature transition is given by 

m{min[Pjq{l3, m)]) 


Oiv(/d) = 


m 


( 6 ) 


N 


Figure [3] shows the pseudo order parameter verses inverse temperature for different 
walk lengths. It can be seen that the order parameter registers a jump from 0 to 1 and 
results in an unmistakable location and estimation of theta temperature. The pseudo order 
parameter is 0 for large value of (3 (small value of temperature). This means the minimum of 
canonical probability distribution occurs for lowest value of contacts, which implies that the 
conformations with low value of contact is least possible. Thus the polymer chain modeled 
by ISAW is in a globule phase, where the number of contacts is large. On the other hand, 
for small value of (3 (high temperature) the pseudo order parameter is 1, implying the high 
contact conhguration least posible. This corresponds to coil phase. 

It has been observed that the inverse theta temperature corresponding to different walk 
lengths fluctuates with N. These fluctuations are observed to be either bound or diminish 
with increasing N (fig. 0]). As compared to earlier study2^, we find that such a fluctuation 
is small on FCC lattice because large nearest neighbour sites. This fluctuation is diminished 
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FIG. 4. Fluctuation of theta temperature for different walk lengths. 

further even for walk length below 100. Thus we limited our analysis for walk length up 
to 100 as mentioned earlier. The inverse temperature corresponding to the theta transition 
estimated using pseudo order parameter is /d© = 0.1385 ± 0.0085. This is in agreement with 
earlier reported result^ of Tg = 7.614 which in terms of inverse temperature is /de = 0.1313. 
Merit of the pseudo order parameter method on FCC lattice is that an accurate estimate of 
theta temperature can be obtained from ISAWs of relatively small walk lengths. In order 
to corroborate our result we also used partition function zeroes method to identify theta 
temperature. 


B. Partition Function Zeros 

In the thermodynamic limit, Yang and Lee proposed a general theory of phase transition 
based on the zeros of grand canonical partition function on complex fugacity plane, which 
they illustrated for Ising model in external held^*^. Fisher extended this approach to the 
zeros of Canonical partition function^. Peter Borrman et ah, showed that it is indeed true for 
systems of hnite size^l. Partition zeros has been used in the study of conformational statistics 
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FIG. 5. Partition function zeros for different walk lengths on the positive quadrant of complex Y 
plane. 

of coil to globule transition of linear homopolymers^‘^'i^‘i2i>2^‘2£. In the thermodynamic limit 
locus of zeros forms a unit circle and crosses the Real axis. The zero which lies on positive 
real axis give the phase transition temperature. But for hnite system they never cross real- 
axis. First zero which is close to positive real-axis in upper half plane becomes important, 
and describes the phase transition temperature. 

The canonical partition function at a particular inverse temperature (3 is dehned in Equa¬ 
tion m which can be recast as 

yy^max 

Cm{P) = E (7) 

m=0 

where y = exp(/d). Partition function is thus a polynomial of order and will have same 
number of roots which are complex. Figure 0 shows roots of the partition function zeros in 
the hrst quadrant of the complex y plane for various walk lengths. 

The hrst zeros does not give us the true value of phase transition temperature because 
of hnite size ehects. So we have to do extrapolation of hrst zeros for longer walk lengths. 
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FIG. 6. Extrapolation of the real part of first zeros of partition function in the limit N ^ oo. 

In order to get a better approximation, we have used slightly bigger walk lengths namely 
100 to 300. The real part of the hrst zeros as a function of 1/N is shown in £g. [6l Let 
Yc be the real part of the hrst zero in the limit N —> oo . The value of /S© = logYc = 
0.1367 ± 0.0049 computed in this method is in agreement with the value obtained using 
pseudo order parameter as well as with earlier reported resultr^. 


IV. CONCLUSION 

Using the notion of atmosphere of polymer growth algorithm we computed the density 
of ISAW on FCC lattice. We estimated the theta temperature of ISAW on FCC lattice by 
using pseudo order parameter method. The theta temperature obtained from pseudo order 
parameter method huctuates around the true value for different walk lengths. Compared 
to other lattices (where the coodination number z is less) we hnd that such a huctuation 
is small for FCC lattice (for which is large). Thus, our pseudo order parameter method 
provides a reliable estimate of theta temperature even for shorter walk lengths. Our result 
showed that the theta temperature of ISAW on FCC lattice obtained from our method 
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agrees very well with the earlier reported result of 7.614. The same value obtained from 
Partition function zeroes provides further support to our theta temperature estimate. 
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